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The  5g  levels  of  atomic  nitrogen 


Edward  S  Chang  and  Hajime  Sakai 


Department  of  Physics  and  Aatronomy,  Univeraity  of  Massachusetts,  Amherst  MA  01003, 
USA 


27  February  1981 


Abstract.  Energy  levels  and  relative  strengths  (or  the  4f-3g  transition  in  N  l  are  calculated  in 
the  /-/  coupling  scheme.  Excellent  agreement  it  found  with  our  interferometric  data  from 
<r-2455  to  tfSoScqi^ in  both  positions  and  intensities.  f 


The  Sg  levels  in  atoms  are  usually  determined  from  measurements  of  the  4f-5g  lines. 
Spectroscopic  studies  were  performed  on  closed-core  atoms  such  as  the  alkali  metals 
(Litzen  1970a)  and  on  helium  (Litzen  1970b).  However,  other  atoms  were  not  studied, 
apparently  because  of  the  expected  complications  in  the  term  analysis  due  to  an  open 
core,  and  in  detecting  low-intensity  atomic  emissions  among  molecular  bands. 

In  this  letter,  ws  show  that  the  Sg  levels  of  N I  (which  has  an  open  core)  are 
accurately  described  by  H  coupling  (i.e.,  in  the  jUCMk  scheme).  This  coupling  ignores 
the  spin  of  the  Rydberg  electron  and  is  an  approximation  to  j-K  coupling  (i.e.,  the 
jKJMj  scheme)  employed,  for  example,  in  11  ui  (Edldn  and  Svensson  197S,  Goldsch¬ 
midt  1981).  The  resulting  simplification  allows  the  analytic  calculation  of  the  relative 
line  strengths  for  the  4f-5g  transition.  This  theoretical  analysis  enables  us  to  identify 
our  glow  discharge  interferometric  data  in  nitrogen  (mixed  with  helium)  from  2455  to 
2S0S  cm-1  as  4f-5g  emission  lines  in  N  i.  These  atomic  lines  lie  in  the  same  spectral 
region  as  the  molecular  W-B  (4, 2)  band  observed  by  Benesch  and  Saum  (1971)  under 
different  discharge  conditions. 

The  general  theory  of  coupling  between  an  open  core  and  Rydberg  electron  with 
high  orbital  angular  momentum  has  been  described  by  Shortley  and  Fried  (1928)  and  by 
Racah  (1942a,  b).  In  application  to  nl  and  nginNn  (Eriksson  1956)  and  to  nf  in  the 
isoelectronic  O  n  (Eriksson  1961),  it  was  noted  that  the  spin  of  the  Rydberg  electron  r 
was  in  effect  decoupled  from  the  rest  of  the  system.  Consequently,  the  system  can  be 
described  in  L-S„  in  H,  or  generally  in  the  intermediate  coupling  scheme  (called  pair 
coupling). 

To  expedite  the  calculation  of  line  strengths,  we  work  entirely  in  the  M  coupling 
scheme.  Since  the  electrostatic  energies  are  usually  given  in  a  coupling  scheme  resulting 
in  total  orbital  angular  momentum  L,  we  apply  the  transformation  to  j-l  coupling  given 
by 


UlLJL,  (LrSe)//*,-(-l),+£'+***K[(2/  + 1)(2/+ 1)]‘/2  { ^  £  ^}.  (1) 


In  equation  (1)  L„  S„  and  /  refer  respectively  to  the  orbital,  spin  and  total  angular 
momenta  of  the  core,  and  /  and  L  are  the  orbital  angular  momenta  at  the  Rydberg 


0022-3700/81/120391  +05S01.50  ©  1981  The  Institute  of  Physics  L391 


TM  U.S.  Qovornmont  li  authorized  to  r •product  and  mi  mu  report. 
Permission  for  further  reproduction  by  others  muit  bo  obtained  from 
Uw  copyrtpbt  owner. 


L392 


Letter  to  the  Editor 


electron  and  the  entire  atom.  AT  is  the  invariant  angular  momentum,  ignoring  s. 
Retaining  only  Fi  of  the  electrostatic  energies  given  by  Racah  (1942a): 


(LS«l£|l5e)-FI(2p.«/)x 


- 1(2/  - 1) 

(2/-1X2/+3) 

-(/+l)(2/+3) 


forL-i+l 
L  -  / 
Z.-/-1 


we  find  the  energy  in  H  coupling  to  be  given  by 

£**-X(2L  +  l)(2/+l)(^  £  {LSAE\LSc) 


(2) 


(3) 


where 


J/i./j./j] 

'/«,/*./•! 

are  the  6- j  symbols.  Using  standard  Racah  algebra,  we  obtain  the  line  strength  as 

SULK,  frm  -  WOK  +  1X2X” + 1)|*  *'  J}*S„  (4) 

where  are  essentially  the  hydrogen  line  strengths,  for  a  Rydberg 

electron  with  I* 3. 

In  application  to  N  I,  we  have  "  S.  - 1  for  the  ground  core  state  (2p)z  *P,  which 
splits  into  /  “0, 1  and  2  levels.  Further  splittings  due  to  electrostatic  interaction  with 
the  7*.  ijerg  electron  are  obtained  from  equation  (3).  The  results  for  the  41  and  the  5g 
configurations  ate  given  in  table  1.  In  addition,  line  strengths  given  by  equation  (4)  are 
shown  for  allowed  transitions  which  obey  A/*0  and  Ait- 0  or  *1.  In  comparing 
relative  strengths  in  H  i,  Bethe  and  Saipeter  (1957)  observed  that  the  most  intense  lines 
are  those  with  An-A/--1.  Our  results  in  table  1  and  die  more  general  results  in 


T*Mt  t.  Theoreticsl  anrjpss  sod  Km  nr«mU»  ot  th«  ♦!«-  5g  trmwrton. 


/ 

K 

ffwc* 

«/4R/3Jf)» 

JC'-l  2  3  4 

5 

0 

3 

0 

0 

1.00 

t 

2 

6 

3 

-7.5 

27.5 

0.71 

0.04 

0.00 

4 

2.3 

-34.5 

0.04 

0.04 

14 

1.22 

2 

1 

-12 

2 

-S3 

0.43  0.12 

0.01 

3 

5.5 

-5.5 

0.40 

0.17 

0.01 

4 

7.3 

32.5 

0.S2 

0.14 

0.01 

S 

-3 

35 

1.10 

0.12 

4 

-24 

1.44 

t  te  mMssT  Me.  41)  -  MM  «w*\ 

Unwa»^FM3a,St)-A0M«?an-1. 

tlemtnstiW. 


Letter  to  the  Editor 


L393 


equation  (4)  enable  us  to  extend  their  rule  to  a  complex  atom  in  j-l  coupling,  i.e.  the 
most  intense  lines  are  those  with  AX  ■  An  *  A/  *  -1.  As  expected  from  properties  of 
the  6-j  symbols,  the  sum  of  the  strengths  for  each  value  of  /  is  (2/  +  1).  The  values  of  the 
electrostatic  interaction  integrals  for  N  t  were  calculated  (Eriksson  1967)  to  be  F2(2p, 
4f)  =  1.861  cm-1  and  F2(2p,  5g) »  0.086  67  cm'1.  In  the  case  of  the  4f  configurations, 
energy  levels  calculated  from  table  1  may  differ  from  the  experimental  values  by  several 
cm-1.  However,  for  the  Sg  configuration  F2  is  substantially  smaller,  so  our  approxima¬ 
tion  can  be  expected  to  be  much  better  as  we  will  see. 

The  infrared  emission  data  analysed  in  the  present  study  were  generated  in  a  ‘Pfund’ 
type  optical  cell  30  m  long,  and  1  m  in  diameter.  At  both  ends  were  placed  two  concave 
mirrors  of  30  m  focal  length.  The  glow  discharge  was  formed  between  the  wall  of  the 
optical  cell  and  a  12  m  long  electrode  placed  at  the  central  region.  An  ac  60  Hz  voltage 
of  700  V  rms  was  applied  to  activate  the  discharge  at  the  rate  of  200  W.  The  source  was 
a  continuous  flow  type;  a  gas  sample  of  a  He-N2  mixture  was  admitted  at  one  end  of  the 
cell  and  pumped  out  at  the  other  end  at  a  constant  rate.  The  pressure  inside  the 


TaM*  2.  Experimental  energy  leveii  of  the  4f  and  the  Sg  configurations  of  N  i. 
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discharge  chamber  was  maintained  at  approximately  0.1  Torr.  The  infrared  radiation 
was  admitted  to  an  infrared  interferometer  placed  at  the  exit  port  of  the  optical  cell 
through  a  KBr  lens.  The  experimental  set-up  was  designed  to  achieve  two  objectives  (i) 
the  generation  of  a  large  number  of  infrared  photons  by  forming  a  large  glow  discharge 
column;  and  (ii)  optimum  spectrometric  efficiency  realised  by  using  the  Fourier 
spectroscopy  technique.  The  interferogram  signal  was  detected  by  an  InSb  detector 
housed  in  a  liquid  N2  Dewar.  The  path  difference  of  the  interferometer  spectrometer 
was  monitored  by  the  interference  fringe  signal  of  the  HeNe  cw  laser  line  6328  A.  The 
spectrum  used  in  the  present  analysis  was  produced  with  a  resolution  of  0.12  cm'1.  The 
accuracy  for  the  present  line  position  measurement  was  considerably  higher  than  that 
obtained  by  the  conventional  spectrometer  measurement,  because  each  spectral  fine  in 
our  measurment  was  directly  referenced  to  the  HeNe  laser  line  (15798.00  cm'1  in 
vacuum). 
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fipn  1.  («)  Our  obsenud  spectrum  between  2455  end  2505  cm"',  (b)  Theoretical 
intoniittaofNiMCsb0wnastidlliMSMlbllo*s.A:/«2,JC«3<-3;B:/>2,JC  -4«-5;C: 
l-l  K »3«-4;  D;  1-2,  K-2-l,  E:  /-l.  K-2«-3;  F:  /- 1.  tC-4-5;  G:  1-0. 
K- 3»4;  H:  /»2,  K-5+6-,  I:  l-l,  lf»3*»4,  and  J:  1-2,  JC «  '•-2.  In  order  to 
dfetineiikh  than  from  the  N I  spectrum,  the  following  He  I  ttnet  are  ihown  as  broken  linet. 
K:4fTF«-5g,Oand4f,F~S|’G.L:4d,D«-Sf,Fand  M:4d,D^5f,F. 


Figure  Ha)  shows  the  emotion  intensity  in  a  nitrogen-helium  glow  discharge.  The 
most  intense  peak  centred  at  er»  2469.58  cm'1  is  due  to  the  4f-5g  line  of  He  i  as  its 
relative  value  varies  with  the  mixture  ratio  of  He  in  several  different  runs.  Note  the 
significant  discrepancy  with  the  previously  measured  value  of  2469.749  cm  ’  (Litzen 
1970b).  We  assign  the  next  largest  peak  at  2489.26  cm'1  as  the  4f-5g  line  (;'  *  2, 
K  *  5«-6)  of  Nt,  which  is  the  most  intense  as  seen  in  table  1.  From  the  experimental 
value  of  the  4f  (/  «  2,  JC  -  5)  level  at  HO  473.37  cm'1  labelled  G(5)  by  Eriksson  and 
Johansson  (1961),  we  determine  the  energy  of  the  5g  (/■  2,  K  -6)  level  to  be 
112962.63  cm'1.  Prom  table  1,  the  rest  of  the  5g  levels  can  be  calculated  by  adopting 
the  experimental  values  for  the  4f  levels.  The  resulting  positions  and  intensities  (<r*  x 
strengths)  are  shown  in  figure  1(6).  The  striking  agreement  between  the  two  sets  serves 
to  identify  the  data  unambiguously  as  the  4f-Sg  transitions  of  N  i.  Simultaneously,  it 
verifies  the  validity  of  j-i  coupling  for  both  the  4f  and  the  3g  configurations.  The 
theoretical  5g  levels  determined  from  table  1  differ  from  our  measured  values  by  only  a 
fete  tenths  of  a  wavenumber. 

fat  conclusion,  the  Sg  levels  of  N  i  have  been  identified  for  the  first  time  and 
measured  to  an  absolute  accuracy  of  0.01  cm'1  by  infrared  interferometry.  The 
emission  line  intensities  to  the  4f  levels  are  in  excellent  agreement  with  theory  using  the 
H  coupling  scheme.  The  general  expression  for  the  line  strength  and  the  propensity 
rate  AJC  --1  are  applicable  to  other  Rydberg  states  in  any  atom  provided  that  /»3. 
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Infrared  emission  spectroscopy  of  glow  discharge  formed  in 
low  pressure  atmospheric  gases 

H.  Sakai,  P.  Hansen,  M.  Esplin,  R.  Johansson,  M.  Peitola,  and  J.  Strong 


Laboratory  study  was  conducted  on  the  molecular  and  atomic  emissions  which  are  considered  detectable 
in  the  IR  radiative  background  of  the  upper  atmosphere.  The  IR  spectra  of  glow-  discharge  emission  formed 
in  air  and  other  atmospheric  constituent  gases  of  0. 1-Torr  pressure  and  a  .'tO-nt  long  column  were  surveyed 
using  the  technique  of  Fourier  spectroscopy.  Several  features  hitherto  unobserved  were  detected  in  our  sur¬ 
vey  study. 


I.  Introduction 

Recently,  in  parallel  with  the  advent  of  IR  technology, 
various  study  efforts  on  IR  emission  revealed  the  atomic 
and  molecular  features  poorly  known  in  the  past.  In 
two  specific  areas  of  the  observational  studies,  one  for 
the  atmospheric  interest  and  another  for  the  astro¬ 
physics  interest,  the  spectroscopic  data  of  the  IR 
emission  are  obtained  for  situations  which  were  con¬ 
sidered  impossible  a  decade  ago.  Most  of  the  emission 
data  currently  obtained  in  field  observations  are  those 
of  a  relatively  low  spectral  resolution,  causing  the  task 
of  making  unambiguous  identification  of  the  spectral 
feature  found  in  the  data  rather  difficult.  A  primary 
objective  in  the  field  observation  study  is  concerned 
with  physical  and  chemical  processes  involved  in  the 
emission  mechanism.  It  cannot  be  achieved  unless  a 
proper  identification  can  be  assigned  to  the  observed 
spectral  feature.  Under  the  circumstances  where  the 
field  data  are  provided  with  a  spectral  resolution  in¬ 
sufficient  for  producing  a  clear  spectral  identification, 
a  laboratory  spectroscopic  study  on  the  candidate 
atomic  and  molecular  IR  emission  provides  worthwhile 
assistance  to  the  analysis  effort.  Our  study  effort  was 
made  in  part  to  fulfill  the  need  called  for  in  such  a  cir¬ 
cumstance.  We  collected  the  data  on  the  IR  emission 
feature  of  atmospheric  species,  atoms,  and  molecules, 
with  a  moderate  spectral  resolution  using  the  technique 
of  Fourier  spectroscopy,  and  studied  them. 
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The  primary  interest  of  our  study  was  focused  on  the 
region  of  the  upper  atmosphere  above  60  km,  where  the 
pressure  remains  below  1  Torr.  The  measurement  ef¬ 
forts  for  the  IR  emission  in  the  upper  atmosphere  were 
summarized  in  a  recent  review  article  by  Baker  et  al. 1 
The  IR  emission  generated  by  the  electrons  of  20-30  eV 
exciting  the  atmospheric  molecules  was  our  main  study 
subject.  Our  basic  experimental  approach  to  the  study 
was  twofold:  (1)  by  forming  a  large  electric  glow  dis-  ; 

charge  column,  the  IR  photons  of  various  atmospheric  | 

species  were  produced  in  a  large  quantity;  and  (2)  by 
using  Fourier  multiplex  spectroscopy,  the  efficiency  of 
collecting  the  generated  IR  photons  was  optimized.2 
By  combining  these  two  factors,  we  were  able  to  enhance 
the  detection  sensitivity  to  such  a  level  that  only  15-roin 
observation  time  is  required  to  cover  the  entire  InSb 
spectral  range  (1800-7800  cm-1)  with  a  spectral  reso¬ 
lution  of  1  cm-1. 

An  atomic  or  a  molecular  transition  falls  into  the  IR 
range  if  the  upper  and  lower  state  of  the  transition  are 
separated  by  an  IR  frequency.  An  electric  dipole 
transition  of  a  homonuclear  diatomic  molecule  takes 
place  from  the  vibrational-rotational  level  of  an  elec¬ 
tronic  state  to  another  level  which  belongs  to  a  different 
electronic  state.  Since  in  general  the  stable  diatomic 
molecules  have  their  first  electronic  excited  state  a  few 
electron  volts  above  the  ground  state,  the  IR  emission 
of  a  homonuclear  diatomic  molecule  is  produced  by  an 
electronic  transition  which  does  not  involve  the  ground 
state.  It  must  take  place  between  two  different  elec¬ 
tronic  excited  states.  The  IR  atomic  emission  occurs  I 

in  a  similar  situation.  It  takes  place  between  two  ex-  t 

cited  states  which  are  separated  by  the  IR  energy. 

These  two  situations  for  the  electronic  transitions  ; 

contrast  to  the  normal  IR  transition,  which  is  found  ’ 

either  in  a  heteronuclear  or  in  a  polyatomic  molecule. 

It  involves  a  vibrational-rotational  transition  within  the  j 

same  electronic  state,  usually  within  the  ground  elec-  L 

tronic  state.  i 
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N.  Experimental  Setup 

A  central  feature  of  our  experimental  setup  is  a  30-m 
long,  1-m  diam  cylinder  used  as  a  container  of  the  glow 
discharge  source,  which  is  formed  between  a  12-m  long 
central  electrode  and  the  external  wall  as  shown  in  Fig. 
1.  An  ac  60-Hz  voltage  of  up  to  1000  V  is  applied  be¬ 
tween  the  electrodes,  activating  the  discharge.  The 
interferometer  accepts  the  IR  radiation  through  a  KBr 
lens  placed  at  the  exit  port  of  the  discharge  source.  The 
interferogram  signal  is  detected  by  an  InSb  detector 
housed  in  a  liquid  nitrogen  Dewar.  The  path  difference 
is  monitored  by  the  interference  fringe  of  the  He-Ne  cw 
laser  line  of  6328  A  (air  wavelength).  The  detector 
output  is  properly  ac  amplified,  synchronously  demo¬ 
dulated,  integrated,  and  finally  converted  to  a  digital 
signal  by  an  analog-to-digital  converter,  which  is 
triggered  by  the  zero  crossing  position  of  the  laser  in¬ 
terference  fringe  reference  signal.  The  digitized  in¬ 
terferogram  signal  is  then  recorded  on  a  mass  storage 
device  (a  floppy  disk)  under  the  control  of  an  LSI/11 
minicomputer.  After  completion  of  the  interferogram 
measurement,  the  interferogram  data  are  postprocessed 
using  our  central  site  large-scale  computer,  CDC 
CYBER  system,  for  the  Fourier  transformation  etc. 

The  optical  cell  which  is  used  to  house  our  glow  dis¬ 
charge  source  is  generally  known  as  the  Pfund  cell.3 
The  discharge  column  formed  between  the  electrodes 
is  seen  thrice  along  the  optical  path,  thus  forming  an 
equivalent  36-m  long  discharge  column.  The  discharge 
region  extends  only  %  of  the  entire  cell  length,  leaving 
%  of  the  cell  a  discharge-free  space.  This  apparent 
disadvantage  created  by  the  discharge-free  region  is  a 
blessing  in  disguise  from  the  experimental  point  of  view, 
since  it  prevents  the  mirror  contamination  caused  by 
the  electron  bombardment. 

An  important  parameter  which  controls  our  glow 
discharge  is  the  excitation  energy  released  to  the  atoms 
and  the  molecules.  A  rough  estimate  of  the  excitation 
energy  can  be  obtained  in  the  following  way.  The 
electric  field  in  the  cell  is  given  as  a  function  of  distance 
r  (measured  from  the  center)  by 


where  V  is  the  electric  potential  applied  between  the 
electrode  of  radius  a  and  the  wall  of  radius  b.  The 
electrons  pick  up  their  energy  under  the  applied  electric 
field  and  lose  it  when  they  collide  with  a  molecule  or  an 
atom.  The  lose  of  energy  occurred  as  the  collision  turns 
into  the  energy  for  exciting  the  collision  partners,  for  our 
case  either  molecule  or  atom.  The  energy  lose  suffered 
at  the  collision  is  usually  resupplied  between  collisions. 
The  mean  free  path  (x )  of  the  electrons,  i.e.,  the  mean 
distance  between  the  successive  collisions,  is  generally 
characterized  by 


<*>  -  !/(«*).  (2) 

where  n  is  the  number  density  of  the  colliding  molecules 
and  <r  is  a  quantity  called  the  collision  cross  section.  If 


Fig.  1.  Schematic  representation  of  the  discharge  column  (not  to 
■cale).  Mirron  Mi  and  Mi  are  focused  on  each  other,  resulting  in 
throe  peases  through  the  discharge  column. 


an  electron  moves  parallel  to  the  field,  it  gains  its  energy 
by  e(x)E,  which  is  calculated  by 


« 


t{x)E 


eVp0 


(3) 


for  the  gas  pressure  p .  In  this  expression,  N0  is  the 
number  density  at  the  standard  pressure  po-  It  is  seen 
that  the  excitation  is  the  highest  in  the  vicinity  of  the 
electrode  and  that  it  reduces  toward  the  outer  wall.  For 
a  typical  example,  we  can  assume  a  *  10-16  cm2,  V  * 
1000  V,  and  n  ■  NoP/po  m  7.07  X  1015  for  p  «•  0.2  Tore. 
The  field  atr  “  30  cm  is  — 13  V/cm,  while  the  mean  free 
path  is  ~1.4  cm.  The  energy  gain  per  collision  in  the 
area  at  r  -  30  cm  would  reach  18  eV  «  145,000  cm-1.  In 
reality,  the  croaa  section  a  varies  dependent  on  the 
electron  energy,  and  the  electrons  may  lose  only  a  part 
of  their  energy  by  the  inelastic  collision  process. 
Nonetheless,  the  value  estimated  above  for  the  excita¬ 
tion  energy  agreed  well  with  the  experimental  data 
obtained.  In  the  glow  discharge  of  oxygen,  for  example, 
observed  are  many  IR  atomic  oxygen  lines  which  are 
produced  in  transition  between  two  highly  excited 
atomic  levels.  By  identifying  the  levels  involved  for  the 
observed  IR  atomic  Ol  lines,  we  estimated  the  excitation 
energy  in  our  glow  discharge  generally  reaching  a  level 
of  20  eV.  Even  though  some  excitations  may  exceed  the 
value  by  a  substantial  degree,  it  is  safe  to  assume  that 
a  majority  of  the  IR  photons  observable  in  our  experi¬ 
ment  are  indeed  generated  in  the  processes  which  re¬ 
quire  an  excitation  energy  of  ~20  eV  or  less. 

One  thing  noticed  is  that  the  excitation  is  very  sen¬ 
sitive  to  the  gas  pressure.  Once  the  gas  pressure  is 
above  0.5  Torr,  the  glow  discharge  which  is  indicative 
of  transitions  between  the  electronic  states  is  confined 
to  the  vicinity  of  the  electrodes,  leaving  a  dark  space 
elsewhere.  For  heteronuclear  diatomic  molecules,  the 
IR  emission  does  not  necessitate  the  electronic  transi¬ 
tion.  An  absence  of  the  visible  glow  discharge  does  not 
eliminate  a  possible  IR  emission.  However,  for  hom- 
onuclear  molecules,  the  visible  glow  discharge  which 
usually  indicates  a  presence  of  the  electronic  transitions 
to  the  ground  state  is  necessary  to  produce  IR  transi¬ 
tion. 

The  glow  discharge  is  run  at  60-Hz  ac  obtained  from 
the  ordinary  60-Hz  power  through  a  step-up  trans- 
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former.  Direct  current  voltages  failed  to  produce  a 
stable  glow  discharge.  The  central  electrode  is  usually 
water-cooled. 

The  spectrum  shown  in  Fig.  2  is  a  typical  result  pro¬ 
duced  by  a  glow  discharge  in  0. 1  Torr  of  air.  Identifi¬ 
able  molecular  features  in  the  spectrum  are  NO,  N2,  CO, 
C02,  N20,  OH,  NH,  and  NO2.  The  atomic  features 
contain  O  and  N.  Each  spectral  feature  of  these  species 
as  well  as  of  others  was  studied  by  observing  the  IR  glow 
discharge  emissions  formed  in  various  gas  samples  in 
addition  to  those  produced  in  air.  Specific  mixtures 
were  selected  to  produce  the  spectral  data  necessary  for 
interpreting  those  features  observable  in  the  air  dis¬ 
charge. 

III.  Emission  Spectra  of  Various  Atmospheric 

Spsciss 

The  spectrum  shown  in  Fig.  2  was  taken  with  a  reso¬ 
lution  of  ~1.0  cm- 1.  A  feature  observable  between  2000 
and  2400  cm-1  consists  of  CO,  N20,  C02,  and  possibly 
others.  The  vibrational  fundamental  (Av  =  1)  of  the 
NO  ground  state  is  observable  in  1800  cm-1,  the  lowest 
frequency  range  covered  by  the  InSb  detector.  A  fea¬ 
ture  which  appears  as  a  group  of  lines  at  2500  cm-1 
consists  of  many  atomic  lines.4  The  transition  fre¬ 
quency  of  2470  cm-1  corresponds  approximately  to  the 
energy  difference  between  n  *  4  and  n  *  5  orbit  of  the 
hydrogenlike  atoms.  Consequently,  all  atoms  con¬ 
tained  in  the  gas  samples  would  emit  IR  lines  in  a  con¬ 
fined  spectral  region  centered  at  2470  cm-1.  We  found 
atomic  lines  of  H,  N,  O,  and  He  in  this  region.  The 
features  observable  between  2500  and  4000  cm-1  are  of 
N2,  NH,  OH,  and  N02-  The  features  above  5000  cm-1 
are  predominantly  those  of  the  N2  first  positive  band 
( R  'ir h~A  32u  ).5  In  addition  to  these  molecular  features, 
various  atomic  lines  of  H,  O,  and  N  spread  over  the  en¬ 
tire  InSb  region. 

The  following  description  provides  a  survey  of  the 
spectral  features  produced  by  various  atmospheric 
species.  As  mentioned  above,  the  glow  discharge 
formed  several  molecular  products  which  were  not 
found  in  the  original  gas  samples.  The  present  paper 
will  not  discuss  the  formation  mechanism  of  these 
products.  All  the  measurements  were  done  in  inte- 
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grating  the  emission  intensity  using  the  normal  lock-in 
technique.  The  formation  processes  of  the  glow  dis¬ 
charge  by-product  would  be  better  studied  by  using  the 
time-resolved  spectroscopic  technique. 

A.  Atomic  Feature 

1.  Hydrogen 

The  atomic  lines  of  hydrogen  were  observed  in  the 
glow  discharge  emission  of  air  whenever  the  sample  was 
moist.  Table  I  lists  the  atomic  H  lines  observed  in  the 
IR.  No  molecular  transitions  of  hydrogen  were  ob¬ 
served  in  the  IR,  even  though  several  known  excited 
electronic  states  of  H2  are  spaced  by  the  IR  frequencies. 
The  atomic  hydrogen  may  have  been  easily  formed 
because  of  a  low  dissociation  energy  at  the  ground  state 
(X'2+).* 

The  glow  discharge  emissions  were  to  a  great  extent 
affected  by  admission  of  hydrogen  into  the  gas  sample. 
There  were  two  notable  effects  observed  in  the  dis¬ 
charge  emission:  suppression  of  the  CO  vibrational 
fundamental  which  would  otherwise  accompany  as  a 
parasitic  emission;  formation  of  two  hydrides,  OH  and 
NH. 

2.  Oxygen 

The  atomic  lines  of  oxygen  were  observed  over  a  wide 
spectral  range.  The  line  position  is  listed  in  Table  II. 
It  is  similar  to  the  hydrogen  case  where  no  molecular 
transitions  were  observed.  The  low  dissociation  energy 
of  the  02  ground  state  (X32^)  contributes  to  the  for¬ 
mation  of  the  atomic  oxygen  in  the  glow  discharge.7 
The  recent  data  compiled  for  the  atomic  IR  lines  by 


Tafcto  L  Obaanad  Atomic  Hydrogen  Unoo  (cm  ’ ) 


An  *  1 

An  ”  2 

An  -  3 

An  -  4 

53.31. 55 

7799.29 

2467.75 

: 1808.25 

2148.79 

or. 

i) 


i  ■  i  »  i  *  i - « - 1  «  r 

2000  4000  6000 

Pig.  3.  Glow  discharge  spectrum  of  Oj. 


Outred4  do  not  contain  the  oxygen  lines  below  3819 
cm-1.  Figure  3  shows  the  glow  discharge  emission 
spectrum  of  the  oxygen  gas. 

Figure  4  is  a  sketch  showing  the  atomic  oxygen  (01) 
energy  levels'*  together  with  the  IR  lines  observed  in  our 
measurement.  The  energy  of  these  levels  is  referenced 
to  the  ground  state  of  atomic  oxygen  (lS)J(2S)2(2p)4 
3Pz,  which  is  located  at  ~40,000  cm-1  above  the  mo¬ 
lecular  O2  ground  state  X3Xg.  As  mentioned  above,  the 


Table  1.  Observed  Atomic  OsyaenUnea 

(cm-1) 

Transition 

7593.7 

4*sS°-3psP 

6289.5 

5d,0-4paP 

5546.9 

4/®P -3d*D 

5479.4 

4/sP-3d»D® 

3918.9 

6f»P-4dsD 

3876.2 

6/a-4 d*D 

3819.9 

8«*S-5p»/» 

3819.9 

8g*G-iPF 

3770.7 

4dtD-4piP 

3617.2 

Ap'P-is'S 

3455.4 

4p*P-4taS 

3226.9 

3d3D-4p3P 

3021.9 

5*»S-4p»P 

2731.0 

5»  39>-*p3P 

2575.7 

5/sF-4d»D° 

2532.7 

bPF-M'D0 

2477.3 

5 g'-G-APF 

2192.1 

4 \p*P-3d*D* 

2154.6 

7 g*G-5f*F 

2150.5 

if*F-5g*G 

TaMeM. 

Observed  Lins  poeWon,  Observed  and  Calcelated  bdenaby  at 

M  TranaMon  between  <*P) •«  and  (*P)4f  ConAguration 

Position 

Relative  intensity 

;  k- 

K" 

(cm-1) 

Observed  Calculated 

0  4 

3 

2484.9 

1.0 

1.0 

1  3 

2 

2480.6 

0.8 

0.7 

4 

3 

2494.1 

1.0 

0.9 

5 

4 

2481.9 

1.3 

1.2 

2  2 

1 

2500.3 

0.5 

0.4 

:t 

2 

2479.0 

0.6 

0.6 

4 

3 

2469.6 

2.7* 

0.8 

5 

4 

2466.7 

1.3 

1.1 

8 

5 

2489.3 

1.4 

1.4 

•  Blended  with  He  1  line. 

estimate  made  on  the  electron  energy  available  fo-  the 
excitation  agrees  with  the  excitation  energy  measured 
from  the  observed  atomic  oxygen  lines. 

3.  Nitrogen 

The  IR  emission  from  the  nitrogen  gas  contains  few 
atomic  nitrogen  lines.  The  data  compiled  by  Outred4 
contain  no  NI  lines  below  5000  cm-1.  The  feature  at 
2470  cm-1  consists  of  many  lines,  as  listed  in  Table  III, 
with  their  transition  identification.9  Figure  5  shows  the 
atomic  NI  emission  at  2470  cm-1. 


-HOOOOcm'1 


-  80000 

3 - 


3 — 

Fig,  4.  Energy  level  diagram  of  atomic  01  indicating  oliserved 
transition*. 


15  January  1952  /  Vo).  21.  No.  2  /  APPLIED  OPTICS  231 


10 


Fig.  5.  (a)  Our  observed  spectrum  between  2466  and  2505  cm'1;  (b) 
theoretical  intensities  of  NI  are  shown  as  full  lines  as  follows:  (a)  j 
•2,K  =  3  —  3;(b);-2,/f-4  —  5;(c); -2.K -3  —  4;(d)>-2,K 
-2-3;(e)y-l.K-2-3;(f)y-l,K-4-5;(g)j-0,lf-3- 
4;<h»-2,/f-5  — 6;<i>y-l,K-3-4;(j>./-2,K-l— 2.  To 
distinguish  them  from  the  NI  spectrum,  the  following  He  I  lines  are 
shown  as  broken  lines:  (k)  4 f'F  —  5 g'G  and  4 pF  —  5 pG\  (I)  4 d‘D 
—  5flF,  (m)  4d30  —  bPF. 


B.  Molecular  Feature 
l.  Nitrogen 

The  N2  infrared  emission  feature  fills  almost  the 
entire  InSb  region  from  2500  to  7800  cm"1.  There  are 
three  distinctive  N2  electronic  transitions  observable 
in  the  IR  (see  Figs.  6  and  7).5  Various  vibrational 
progressions  of  the  first  positive  band,  are 

observable  over  a  wide  spectral  range.  The  emission 
spectrum  of  the  N2/He  glow  discharge  is  shown  in  Fig. 
8  to  illustrate  the  band  formations.  The  band  centers 
of  B37r„-A3Z,J,10  which  are  indicated  in  the  Figure,  are 
listed  in  Table  IV.  In  contrast  to  the  first  positive  band, 
the  other  two  electronic  transitions,  the  Wu-Benesch 
band  W  'Au-B3jrH  (Ref.  11)  and  the  McFarlane  band 
t^),Au-a,lr<,,2  do  not  form  a  well-recognizable  vibra¬ 
tional  progression.  The  only  observable  transition  for 
the  IV:,AU-BV„  band  is  formed  at  the  3000-cm-1  re¬ 
gion.  The  w'Au-o  'Z*  band  is  only  observable  in  the 
2750-cm-1  region  where  the  Wu-Benesch  (3-1)  band 
should  be  found.  In  the  2470-cm-1  region  where  the 
Wu-Benesch  band  should  be  found,  the  emission  fea¬ 
ture  consisting  of  a  group  of  the  atomic  nitrogen  lines 
is  observed. 

The  excitation  of  the  first  positive  band  B3*>-A3Zj 
is  found  to  be  strongly  influenced  by  the  discharge 
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Fig.  7.  Spectrum  of  the  ut'A-a'r  (0-0)  transition. 


Frequency  cnH 

Fig  8.  Glow  discharge  spectrum  of  air  mixed  with  H|.  Mixture  ratio:  air.  3  and  H*  1. 
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condition;  its  intensity  varies  over  a  large  range  from  an 
undetectably  weak  level  to  an  extremely  strong  level. 
The  other  two  hands  which  are  seen  in  the  low  fre¬ 
quency  region  remained  consistently  strong.  The 
characteristics  of  both  hands,  one  at  2740  cm-1  and 
another  at  :X)00  cm-1,  remain  consistent;  a  singlet 
structure  was  observed  for  the  27 40 -cm'1  band  shown 
in  Fig.  7  and  a  triplet  structure  of  the  3000-cm”1  band. 
The  2740-cm-1  band  was  identified  as  the  u< 1  A*-o  lxf 


transition.  We  were  able  to  determine  the  spectro¬ 
scopic  parameters  of  this  band  and  found  them  in  good 
agreement  with  those  determined  from  the  UV  spectral 
data. 

2.  OH 

It  was  found  that  enhancement  of  the  hydride  emis¬ 
sion,  OH  and  NH,  accompanied  an  increase  of  the 
moisture  in  the  sample  air.  The  spectrum  shown  in  Fig 
8  is  the  glow  discharge  emission  formed  in  the  air  mixed 
with  hydrogen.  Our  experiments  indicate  that  these 
hydrides  are  formed  with  the  presence  of  hydrogen, 
either  atomic  or  molecular.  The  presence  of  HgO  does 
not  seem  to  be  required  for  the  formation  of  the  hy¬ 
drides.  The  infrared  OH  band  was  thoroughly  studied 
by  Maillard  et  al.n  Our  measurement  of  the  OH  line 
position  agrees  well  with  them. 

3.  NH 

A  unique  feature  discovered  in  our  study  is  the  vi¬ 
brational  fundamental  transition  of  the  NH  free  radical 
in  its  electronic  ground  state  (3Z~).  The  NH  lines  were 
observed  in  the  discharge  formed  in  air.  It  was  gener¬ 
ally  considered  that  the  decomposition  process,  NHS  -* 
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Fig.  9.  Infrared  spectrum  of  NH  (s2);  spectral  resolution  0.12  cm*1. 


NH2  -*  NH,  is  the  feasible  generation  mechanism  for 
this  free  radical.  However,  the  spectra  taken  with  the 
air/H  j  mixture  showed  not  only  a  predicted  increase  in 
the  OH  band  but  also  a  striking  enlargement  in  the  NH 
bands.  A  presence  of  NH  in  the  glow  discharge  emis¬ 
sion  found  in  the  air/H2  or  N2/H2  mixture  indicates  that 
the  recombination  process  is  undoubtedly  responsible 
for  NH  formation  in  these  gas  samples. 

The  vibrational-rotational  transition  of  NH  exhibits 
a  more  complex  structure  than  that  of  the  singlet.  The 
NH  spectrum  in  Fig.  9  was  taken  with  a  resolution  of  0.1 
cm~l,  which  was  adequate  for  spectral  analysis.14 
Table  IV  lists  the  spectroscopic  constants,  which  are 
more  accurate  than  those  previously  measured  using  the 
UV  data.15 

4.  NO 

The  NO  band  was  seen  in  the  air  discharge  shown  in 
Fig.  2.  A  complementary  relation  was  observed  be¬ 
tween  the  intensity  of  the  NO  and  the  NH  bands.  With 
introduction  of  extra  hydrogen  to  the  discharge  gas,  a 
weakening  of  the  NO  band  was  observed  with  an  en¬ 
hancement  of  the  NH  band.16 

IV.  Molecular  Emission  In  the  2000-2400-em-1 
Region 

The  feature  observed  between  2000  and  2400  cm-1 
is  composed  of  several  bands,  the  CO  fundamental,  the 
N2O  Ai»n  *  1  transition,17  the  CO2  Av;i  ■  1  transition, 
and  possibly  others.  The  intensity  pattern  was  found 
critically  influenced  by  the  discharge  condition.  The 
most  persistent  feature  is  the  03  transition  of  NgO.  The 
excitation  of  the  CO2  band  is  affected  by  the  presence 
of  N2.  In  our  glow  discharge  condition,  an  exact  nature 
of  the  correlation  between  the  CO2  and  the  N2  excitation 
is  unclear.  A  detailed  study  of  the  emission  in  this  re¬ 
gion  requires  a  spectral  resolution  better  than  0.1  cm-1, 
which  is  the  limit  of  our  experiment. 
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A  NEW  ROTATION-VIBRATION  HAMILTONIAN 
L.S.  Mayants 

Abstract 

The  conventional,  rather  tangled  rotation-vibration  Hamiltonian, 
whose  coefficients  are  expressed  in  terms  of  normal  coordinates,  is  shown 
to  be  appropriate  only  for  small  enough  atomic  displacements.  A  nev, 
exact  rotation-vibration  Hamiltonian  is  obtained  in  an  explicit  form, 
whose  coefficients  can  be  expressed  in  terms  of  Cartesian  components  of 
the  atomic  displacements.  The  method  used,  for  the  construction  of  this 
Hamiltonian  is:  (l)  calculation  of  the  matrix  of  derivatives  of  the  set 
of  atomic  Cartesian  coordinates  in  the  initial  rest  coordinate  system  vith 
respect  to  the  set  of  coordinates  describing  separately  translational, 
rotational,  and  internal  motions  of  a  molecule;  (2)  construction  of  the 
inverse  matrix;  (3)  calculation  of  the  kinetic  part  K*  of  the  rotation- 
vibration  Hamiltonian;  (h)  similarity  transformation  of  K'  by  making  use 
of  the  Jacobian  of  the  initial  coordinates  with  respect  to  the  new  ones, 
which  cancels  the  terms  linear  in  momenta  canonically  conjugate  to  the 
internal  coordinates;  (5)  obtaining  the  final  form  of  the  Hamiltonian. 

The  new  Hamiltonian  allows  one  to  represent  it  easily  in  the  form  of  the 
corresponding  differential  operator,  which  opens  the  door  also  to  the 
utilization  of  methods  of  approximate  computation  of  rotation-vibration 
spectra  other  than  the  perturbation  one.  It  is  proven,  by  the  way,  that 
the  supposition  of  the  existence  of  an  equilibrium  configuration  of  a 
molecule  is  quite  enough  for  constructing  the  Hamiltonian,  Eckart's 
conditions  turning  out  to  be  not  only  approximate,  but  also  unnecessary. 


The  new  rotation-vibration  Hamiltonian  can  prove  to  be  of  special  importance 
for  computation  of  molecular  apectra  in  the  case  of  large  rotation  or 
vibration  quantum  numbers,  in  particular  for  computation  of  spectra  of 
molecules  near  their  dissociative  limits. 

1.  Introduction 

A  rotation-vibration  Hamiltonian  was  introduced  for  the  first  time  by 

1  2 
Wilson  and  Howard,  and  modified  somewhat  later  by  Darling  and  Dennison. 

Its  general  form  was  considerably  simplified  by  Watson. ^  Nevertheless,  the 
conventional  Hamiltonian  still  remains  rather  tangled.  The  fact  that  its 
coefficients  are  expressed  in  terms  of  normal  coordinates,  the  atomic  dis¬ 
placements  being  supposed  to  be  linear  combinations  of  those,  makes  the 
Hamiltonian  approximate,  because  in  the  general  case  the  supposition 
mentioned  cannot  hold.  Besides,  the  same  fact  leads  to  certain  difficulties 
in  the  practical  utilization  of  the  Hamiltonian,  in  particular  Bince  its 
coefficients  are  implicit  functions  of  the  force  coefficients  and  should  be 
changed  in  the  process  of  computation  of  spectra,  if  the  potential  energy 
is  improved. 

The  methods  of  the  numerical  computation  of  the  rotation-vibration 
spectra  based  on  the  conventional  Hamiltonian  can  apparently  be  only  various 
versions  of  the  perturbation  method.  These  versions  usually  utilize  up  to 
the  sixth  order  approximations  of  that  method.  In  the  meantime,  the 
resolution  of  rotation-vibration  spectra  has  been  greatly  increased  due  to 
the  development  of  technique  of  infrared  molecular  spectroscopy.  And  even 
though  the  conventional  methods  of  computation  proved  apparently  to  work 
sufficiently  good  in  many  cases,  they  are  not  aeeurate  enough 


for  the  computation  of  the  spectra  in  the  case  of  rather  large  rotation  or 
vibration  quantum  numbers. 

Thus  the  accuracy  of  the  computation  of  the  rotation-vibration  spectra 
turned  out  to  fall  behind  the  accuracy  of  the  contemporary  experimental 
data.  This  fact  requires  increasing  the  accuracy  of  the  computation  of  the 
spectra.  But  the  higher  the  accuracy  required,  the  higher  should  be  the 
order  of  the  perturbation  method  approximation  to  be  employed.  And  this 
leads,  in  turn,  to  more  and  more  complicated  computations.  Besides, 
appearance  of  even  an  approximate  occasional  degeneracy  in  the  energy  levels 
causes  the  necessity  of  utilizing  some  special  computational  methods.  And 
all  the  difficulties  cannot  be  overcome  when  using  the  conven¬ 
tional  Hamiltonian,  the  more  so  because  this  Hamiltonian  is  not  quite 
accurate  by  itself. 

Under  the  circumstances,  it  seems  interesting  to  try  to  obtain  a 
rotation-vibration  Hamiltonian  in  explicit  form  and  expressed  in  terms  of 
Cartesian  displacements  of  atoms,  which  would  make  it  generally  valid.  The 
aim  of  this  article  is  to  shov  how  this  task  was  accomplished. 

In  order  to  obtain  the  Hamiltonian  needed  it  was  necessary  to  give  up 
the  well  known  common  way  of  obtaining  a  Hamiltonian  in  coordinates  other 
than  Cartesian  ones,  and  to  try  a  different  method.  The  method  used  is  the 
direct  representation  of  the  Hamiltonian  originally  written  down  in 
Cartesian  coordinates,  in  the  coordinates  describing  the  translational, 
rotational,  and  internal  motions  of  a  molecule,  separately. 

In  any  coordinates,  the  Hamiltonian  has  the  form 

H  ■  K  ♦  V  (1) 

where  K  is  the  kinetic  and  V  is  the  potential  part  of  the  Hamiltonian. 


Certainly,  only  K  should  be  considered,  since  V  does  not  need  any  trans¬ 
formation. 


In  the  initial  rest  coordinate  system  (laboratory  system)  the 
differential  operator  K  has  the  fora 


K  ®  \  Z  e.(p^,  p^h. 
i-1 


Here 
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r^  ■  (Xj.yj.Zj)  is  the  radius  vector  of  ith  atom;  Oj  *  («»«x,y,s)  is 

the  oth  Cartesian  coordinate  of  ith  atom;  e.^mT1  i»  the  inverse  mass  of 

ith  atom;  and  n  is  the  nuaber  of  atoms  in  a  molecule. 

We  choose  the  nev  coordinates  as  follows:  X°  •  R°,  Y°  *  R°,  and 

*  * 

Z°  •  R®  are  the  Cartesian  coordinates  of  the  center  of  mass  of  the  molecule 
z 

in  the  initial  rest  coordinate  system;  4  ■  0  •  Qg,  and  X  *  $1^  are 

Euler's  angles  determining  the  spatial  orientation  of  the  molecule  equi¬ 
librium  configuration  and,  hence,  the  spatial  orientation  of  the  rotating 
coordinate  system  rigidly  attached  to  it;  q^  (J*l ,2, . ,.^l)  are  some  indepen¬ 
dent  internal  coordinates  describing  a  deformation  of  the  molecule  relative 
to  its  equilibrium  configuration,  and  H  ■  3n-6  is  the  number  of  these 
coordinates.^ 

In  the  chosen  nev  coordinates,  Eq.  (3)  yields 

±  .  ~  J  - 3 
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By  introducing 


we  get  from  Eqs.  (3)  and  (M 


(5) 


p(1>  ■  *  '*iVj  <« 

(From  now  on  we  will  widely  use  matrix  calculus.  We  will  use  the  same 

notation  a  *  (a^,a,),. . .  ,a^)  for  both  the  vector  ji  vith  components  and 

the  column  matrix  a  with  elements  a.  (j*l,2, . . . ,t) .  A  row  matrix  b  with 

“* 

elements  J>j  ( j*l,2, . . . ,k)  will  be  written  as  £  ■  |  |bjbg. .  .bjJ  | .  cTwill 

denote  the  matrix  transposed  to  £.  The  saw  notation  will  apply  to  block 

matrices.  En  will  denote  an  n  by  n  unit  matrix. ) 

It  only  remains  now  to  find  all  the  VjR®,  ,  and  ,  and  to 

substitute  them  in  Eq.  (6),  and  then  to  replace  the  p^  in  Eq.  (2)  by 

making  use  of  the  results  obtained.  It  is  very  simple  to  calculate  V.R°. 

I  a 

Indeed,  by  definition 


f°  -  M"*1  E  m.r'1'  (7) 

i  1 

where  5°  is  the  radius  vector  of  the  center  of  mass  of  the  molecule  and 


M  ■  Z  m .  is  the  total  mass  of  it.  Therefore 
i 


V1  8°  «  M"1  EyBj  (8) 

where  E^  is  a  3  by  3  unit  matrix. 

But  it  is  not  so  simple  to  find  VjtVj  (j-1,2,3)  and  (j-1,2,. . .  ,N) , 

because  it  is  impossible  to  represent  0,  0,  X,  and  all  the  q.  as  explicit 
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functions  of  the  xt,  y^  and  (i»l,2,. . . ,n).  To  find  these  quantities, 
we  calculate,  first,  the  3n  by  3n  matrix 
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(a,  0  *  x,y,z).  Then  we  construct  the  inverse  a&trix: 
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from  which  we  can  obviously  get  all  the  quantities  we  need,  since 
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The  purely  analytical  calculation  of  matrix  D  is  given  in  f2.  The 
construction  of  matrix  D_1,  by  making  use  of  some  important  relations 
concerning  the  theory  of  small  molecular  vibrations,  is  presented  in  113-5* 
In  Sf6-8,  the  kinetic  part,  K' ,  of  the  rotation-vibration  Hasdltonian  is 
constructed  in  terms  of  angular  and  vibrational  angular  momenta  operators. 
Since  K*  contains  terms  linear  in  Pj ,  which  do  not  disappear  in  the  zero 
approximation,  an  appropriate  similarity  transformation  of  that  canceling’ 
those  terms  is  performed  and  the  final  rotation-vibration  Hamiltonian  is 
calculated  in  119-10.  In  the  calculation,  general  formulas  for  commutators 
of  functions  of  more  than  one  pair  of  canonically  conjugate  operators  are 
used. 
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2.  Calculation  of  Matrix  D 

The  transformation  of  any  vector  a  given  in  the  rest  coordinate  system 
into  the  corresponding  vector  aj  in  the  rotating  coordinate  system  is  of 


the  form 


where  is  the  orthogonal  matrix^ 
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Since  A^-1  =»  the  radius  vectors  r^  and  r^  of  ith  atom  in 


(15) 


o  o 

the  rest  and  rotating  coordinate  systems,  respectively,  are  related  by 

j(i)  „  Jo  +  J«(n)^( i)  ’  _ 

o 

The  3n  Cartesian  components  fia^i-l.a,. . .  ,n;  a*x,y,z),  in  the  rotating 
coordinate  system,  of  the  atomic  displacement  vectors 

sf<l> . 

are  generally  nonlinear  functions  of  any  internal  coordinates  q^  [5^^  is 
the  equilibrium  radius  vector  of  ith  atom  in  the  rotating  coordinate  system]. 
In  virtue  of  Eq.  (7),  the  equation 

/i\i 
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is  valid  for  any  arbitrary  fir^. 


At  the  same  time,  the  veil  known 


equation 
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i 3  valid  only  for  small  enough  fir 
Section  U. 

From  Eq.  (15)  ve  get 
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This  statement  will  be  proven  in 
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Simple  calculations  yield 
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where  the  matrices 
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are  closely  related  to  infinitesimal  rotations  around  the  axes  X' , 

O 

Z'  of  the  rotating  coordinate  system,  respectively. 

Zq.  (20),  in  view  of  Eq.  (22),  takes  on  the  form 
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The  upper  parenthetic  index  "r"  indicates  that  A 


(r)« 


(26) 


is  related  to  rota¬ 


tions;  the  columns  of  this  matrix  being  the  "modes"  of  small  rotations  of 

.(O’ 


around  the  respective  axes  X' ,  Y' ,  and  Z'  of  the  rotating 
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the  vector  r 
coordinate  system 

In  Eq.  (2l),  we  can  make  use  of  the  relation 
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where  A  is  a  one  by  three  column  matrix  whose  elements  are  determined  by 
1 J 

the  choice  of  internal  coordinates,  the  geometry  of  the  equilibrium  con¬ 
figuration  of  the  molecule,  and  masses  of  its  atoms.  We  obtain  thereby 
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We  can  rewrite  A.v  in  Eq.  (25),  in  view  of  Eq.  (l6),  as 
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The  columns  of  the  composite  matrix 
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represent  the  modes  of  small  rot at  ions  of  the  equilibrium  configuration  of  the 

o 

molecule  around  the  respective  axes  in  the  rotating  coordinate  system. 

We  will  also  utilize  the  composite  matrix 

a  *  It  hi  |  ■  (ai . \ . v- 

We  can  now  rewrite  matrix  D  as  given  in  Eq.  (9)  in  the  form 
D-  He’*!  »(fl,A(r,,c(8,j  X(q)a||. 
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the  3n  by  3n  diagonal  matrix  of  the  form 
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and  A'  is  obtained  from  Ax  ;  by  substitution  of  r'  '  for  R'  thus,  in 


view  of  Eq.  (29), 
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where  6A  is  composed  of  5A^  in  the  same  way  as  A  of  A. 


Let  us  note,  by  the  way,  that  the  columns  of  the  matrix  E'  represent  the 
modes  of  the  translations  of  the  molecules  as  a  single  whole  in  the  initial 
rest  coordinate  system. 


,-l 


3.  Construction  of  the  Matrix  D 

To  construct  the  matrix  D-1  as  given  in  Eq.  (11),  we  will  make  use  of 
the  obvious  relation 
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From  Eq.  (8)  we  immediately  obtain 
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(*iO) 


fa0  ■  M-ifeyn 


where 


“"||“iaij|l  (U1) 

and  it  remains  only  to  find  and  ^q.  In  order  to  perform  this  task,  we  will 
use  some  relations  concerning  the  theory  of  small  molecular  vibrations. 


U.  Some  Necessary  Relations 

Let  fir  *  (fir^ . .  ,fir^  ^  , . . .  ,fir^ )  be  the  column  matrix  representing 


the  totality  of  infinitesimal  changes  in  Cartesian  coordinates  of  an  n  atomic 

molecule  in  the  initial  rest  coordinate  system;  fir^  *  (fix^.fiy^fiz^)  being 

the  one  by  three  column  matrix  related  to  ith  atom.  Let,  further,  q^  *  (q^, 
( r ) 

q  ,q)  be  the  column  matrix  representing  the  complete  set  of  coordinates 

describing  separately  translations,  q^\  rotations,  q^ ,  and  internal 

(c) 

motions,  q,  of  the  molecule.  Then  the  relation  of  fir  to  q  can  be  repre¬ 
sented  in  the  two  following  forms 


(c)  (c) 
fir  ■  A  q 


q(c)  a  B(c)  6r> 


the  condition 


B(c)A(c) 


J3n 


(k2> 

(U3) 

(MO 


(c) 


having  to  be  satisfied.  In  accordance  with  the  form  of  q  ,  we  rewrite 
A(c)  and  B(c)  as 

A<c>-  |U(t>  !  A<->  |a||  (1,5) 

B^  *  (B^.B^.S).  (U6) 

Equations  ( 1*2)  and  (U3),  in  view  of  Eqs.  (U5)  and  (U6),  respectively, 

yield 
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+  Aq 


( *»7 ) 


(t)q(t)  +  A(r)q(r)  +  Aq 

(!*7) 

q(t)  -  B«>6r 

(J*8a) 

q(r)  -  B(r)dr 

(U8b) 

q  *  Bdr. 

(U8c) 

It  follows  from  Eq,  (hk )  that 


B«>A(t>  « 

E3 

(»»9a) 

B(r)A(r)  - 

E3 

(»*9b) 

BA-E,, 

(*»9c) 

B(t)A(r)  . 

0 

( i»9d ) 

B(t)A  -  0 

(l»9e) 

- 

0 

(l*9f) 

B*r)A  *  0 

(*»9g) 

BA^  *  0 

(l*9h) 

BA^r)  =  0 

(»»9i) 

From  Eq.  (**7)  we  get 


36r 


(t).  _36r 


A(r);~-A. 

3q 


(50) 


By  comparing  this  with  Eqs.  (19) ,  (2k),  and  (28),  respectively,  we  find,  in 


view  of  the  fact  that  for  the  case  under  consideration 


*(fl) 


J3’ 


A{t)  *  E^  (51) 

while  matrices  and  A  are  those  given  in  Eqs.  (32)  and  (33)»  respectively. 

The  matrix  B  is  the  well  known  matrix  used  in  calculations  of  molecular 
vibrations.  Its  elements  can  be  easily  found  and  depend  only  on  the  choice 
of  internal  coordinates  and  the  geometry  of  the  molecule. 1-0-12  To  find  of 
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what  a  fiirm  me  muLri  l/1'^  and  B^r^,  w«  will  make  uae  of  the  equation 


T(C)  •  *<' WC> 


(52) 


which  connects  the  kinetic  energy  matrix  m  [see  Eq.  (kl)]  in  the  fir  coordi¬ 
nates  with  the  total  kinetic  energy  matrix  T^c ^  in  the  q^  coordinates.  In 

the  case  under  consideration,  translations,  rotations,  and  internal  motions 

(c) 

of  a  molecule  are  fully  separated,  and  we  can  represent  T  ,  in  accordance 

(c) 

with  the  form  of  q  ,  as 


.<«> 


T(t)  0 


0 

0 


0 


(53) 


where  ,  T^r\  and  T  are  the  kinetic  energy  matrices  for  the  translations, 

rotations,  and  internal  motions,  respectively. 

From  Eq,  (52)  we  get.  In  view  of  Eq.  (kk), 

B(c)  .  t<c)-1*(c)b 

Equation  (5k)  yields,  in  view  of  Eqs.  (k5)  and  (k6), 

«(*> . 


where 


,(t) 


B<r)  »  0(r¥r)» 

B  ■  j  |  ...  Bj^ ...  |  J  »  G 2a 


G""  «  G(r)  .  ^(r)-l,  Q  .  j-1^  (56) 

Comparing  Eqs.  (52)  and  (53),  we  find,  in  view  of  Eq.  (51),  and  Eqs.  (32) 
and  (30),  respectively, 

■  fe'mE'  ■  EmtE3  •  ME^  (57) 


(5k) 

(55a) 

(55b) 

(55c) 


and 
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p(r)  -  x(rWr) 


I  - 


where 


I 

aa 


I 

V 

-1 

-I 

XX 

xy 

xz 

-I 

1 

-I 

xy 

yy 

yz 

-i 

-i 

I 

xz 

yz 

zz 

+  I  («4b^y); 

Y 


x., '  jvi11"!1’  <**>• 


(58) 


(58e) 


Hence 

-  M*1^  (59a) 

G*r)  -  I-1  (59b) 

(I  Is  known  to  be  the  moment  of  inertia  matrix.  G  la  the  veil  known  matrix 
used  in  calculations  of  molecular  vibrations ,  whose  elements  for  a  number  of 
internal  coordinates  are  tabulated  in  some  works.10"12) 


Thus  we  finally  find 


»  M-1^m 

(60a) 

B(r)  . 

(60b) 

13 

We  will  also  need  the  equation 

<v 

A  *  eBT 

(61) 

following  ionediately  from  Eq.  (55c)  ,  where 

”*‘l-  ll'i'nll 

(62) 

and  the  equation 

B.S-0-  ||gj,|| 

(63) 

which  results  from  Eqs.  (6l)  and  (49c). 
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We  can  now  easily  see  that  Eq.  (49 f)  Is  a  corollary  of  Eq.  (49d),  and 
Eq.  (49h)  is  that  of  Eq.  (49e);  both  pairs  of  equations  reflecting  the  law  of 
conservation  of  the  center  of  mass  of  the  molecule.  Eq.  (491)  is  the 
corollary  of  Eq.  (49g)  and  both  can  be  rewritten,  in  view  of  Eqs.  (32)  and 
(30),  as  a  set  of  equations 

^bla)ReU  "  bw>Hi1>)  *  0  (a’B"x’y’z)  (6U a) 

or 

lmi(^l)ajj)  -  R^aj^)  -  0  (o,B«x,y,x)  (64b) 

where  b^  is  an  element  of  the  matrix  B  and  a^  is  the  corresponding 
element  of  the  matrix  A;  the  upper  parenthetic  index  i  numbering  the  atoms, 
the  lover  index  J  numbering  the  internal  coordinates,  and  the  index  o 
referring  to  the  «th  Cartesian  coordinate  of  an  atom. 

Equation  (l8)  results  from  Eq.  (64b),  provided  Eq.  (47)  is  valid>  But 
Eq.  (47)  cannot  generally  hold  for  any  4r^ — it  should  be  valid,  generally 
speaking,  only  for  small  enough  4r^  that  can  be  represented  as  linear 
functions  of  any  internal  coordinates  q^.  And  it  is  precisely  this  fact 
which  proves  the  statement  following  after  Eq.  (18). 

Equation  (64)  should,  in  turn,  be  valid  for  any  arbitrary  choice  of 
Internal  coordinates  q^,  which  is  possible  only  because  the  matrices  B  and 
B' ,  respectively,  A  and  A'  for  the  corresponding  different  sets  of  coordi¬ 
nates  q  and  q'  are  related  by  a  linear  transformation.  Hence,  the  fulfill¬ 
ment  of  Eq.  (64a)  for  all  the  matrices  B  follows  from  their  fulfillment  for 
any  of  them.  Matrix  B  has  the  simplest  form  when  small  changes  in  the 
distances  between  pairs  of  atoms  are  chosen  as  internal  coordinates.  In 
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■  vr> 


this  case  the  internal  coordinates  q, . ,  corresponding  to  the  pair  of  ith  and 
i'th  atoms  has  nothing  to  do  with  all  the  other  atoms,  and,  therefore,  all 
the  elements  of  the  respective  rov  of  matrix  B  are  zeros,  except  for  those 
referring  to  jLth  and  ^Jth  atoms,  which  are**1  (up  to  a  sign)  (R^^  -  R^*  ^)l~*, 
and  (R^*  ^  -  R^^l^*,,  respectively,  where  1^,  is  the  equilibrium  distance 
between  ith  and  i'th  atoms.  Substituting  this  in  Eq.  (6Ua),  we  get 

9mm 

identically 


(R^-R^’^R^-fR^-R^'W^+fR^'W^R^-fR^-R^W1'^  (65) 

a  a  p  p  p  a  a  a  p  p  pa 

Thus,  Eq.  (6Ub)  for  any  choice  of  internal  coordinates  follovs  directly 
from  the  identity  represented  by  Eq.  (65)*  But  this  means  that,  in 
distinction  to  Eqs.  (l*9d)  and  (49e) ,  Eq.  (U9g)  does  not  reflect  any  physical 
law.  On  the  contrary,  Eq.  (l8),  which  can  be  interpreted  as  the  equality  of 
the  angular  momentum  corresponding  to  small  deformations  of  a  molecule  to 
zero,  results  in  the  last  analysis  from  the  mathematical  identity  (65) .  And 
this  fact  eliminates  all  the  difficulties  related  to  the  problem  of  choosing 
the  rotating  coordinate  system.*'* 


5.  Construction  of  Matrix  D~*  (continued) 

It  is  obvious  that  $R°  as  given  in  Eq.  (40)  coincides  with  B^  as 
given  in  Eq.  (60a).  Therefore,  Eqs.  (39a),  (39d),  and  (39e),  in  view  of 
Eqs.  (U9a),  (b9d),  and  (b9e),  are  satisfied  (the  commutativity  of  with 
XiQ)  is  to  be  taken  into  consideration).  In  order  that  matrix  satisfy 
Eqs.  (39f)  «nd  (39g)  it  should,  in  view  of  Eqs.  (5l)t  (*»9f),  (l»9g),  and 
(60b),  have  the  form 

■  c'Jt*rWn*  (66) 
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where  C'  is  a  3  by  3  matrix  determined  by  the  condition 

c’*(rWr),c(Q)  .  ^ 

which  follows  from  Jq.  (39b).  Simple  calculations  give 

•  »  » 


XX 

S» 

_I1X 

• 

1 

1 

xy 

1yy 

-I 

zy 

• 

• 

t 

XI 

"V* 

I 

zz 

where 


(a,  0,  y  is  any  of  x,  y,  z).  Therefore 


C*  ■  C^'o^* 


where 


sink 

sine 

cosy 

sln6 

0 

-  c(Q)-1  - 

-  cosx 

sinx 

0 

(71a) 

sinxctge 

cosxctgO 

-1 

I 

0(r)' 

■  II 

1  ■ I''1 

i 

(71b)  ! 

1 

and  fa  takes  on  the  final  form 

fa  ■  c*0*Vr*'tfrW°*.  (1 

If  we  choose  fa  ■  BA^°\  then  Bqs.  (39c)  and  (39h)  would  be  satisfied, 
because  of  Bqs.  (U9c)  and  (49h) .  But  Eq.  (391) will  not  be  fulfilled,  since, 
in  view  of  Eqs.  (37)  and  (U9i), 

BA(r)'c(°)'  ■  B«A*rMn*'  +  9.  (1 
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In  order  not  to  disturb  the  fulfillment  of  Eqs.  (39c)  and  (39h),  and  to  ensure 
the  satisfiability  of  Eq.  (39i)»  we  should  choose  in  the  form 

*4  -  <B  -  xG(r)l{r)m)A(ft) 

where 

*  "  (*1—- •v— V  *  lh«ll 

is  an  H  by  3  matrix  determined  by  the  obvious  conditions 

BA(r)'  -  KG{r),*(r)aA(r)\ 

This  gives  for  x,  in  view  of  Eqs.  (68)  and  (71b), 

x  -  ■  MA^r^.  (77) 

Hence,  the  final  expression  for  $q  is  as  given  in  Eq.  (7*0  • 


(74) 

(75) 

(76) 


6.  Construction  of  the  Operator  K  in  the  Hew  Coordinates 
From  Eqs.  (Uo),  (72),  and  (74)  we  get,  respectively, 

■  M_1miE3 

vs .  x<nv»Sr)s(r)’*(°)’ 

1  011 

This  gives,  in  view  of  Eq.  (6), 

P(i)  -  M-^PO  +  ^0,m1Ajr,8(r)'(2f(Q),H  -  fe)  ♦  ^0)SiP 

where 

P°  »  a“  <WV‘  P  " 

Upon  introducing  Eq.  (8l)  into  Eq.  (2),  we  will  generally  get  terms  both 
quadratic  and  linear  in  P°,  n,  and  P.  In  calculation  we  will  use  the  matrix 
representation  of  the  scalar  product  of  vectors,  namely, 

(a,b)  ■  alb.  (83) 


(78) 

(T9) 

(80) 

(81) 

(82) 
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A^Jpd)  for  thus  obtaining 


m'  •  Er^'o  A*nV^. 


a  o 


(89) 


Since  r  1  °a»  in  viev  of  E<1*  (25),  is  ath  row  of  matrix  ,  we  get  from 

Eq.  (89) 

N  *  lX(r),A  (n)p(1).  (90) 

1  i  o 

This  gives,  because  of  Eqs.  (39i),  (8l),  and  (68), 

#4(r),rt(r)i(rl,«(fl,,n.J(a),n.  (9i) 


Besides,  we  can  Introduce,  in  view  of  Eqs.  (55c)  and  (7T), 

V  *  «  X^’mAGP  ■  <P.  (92) 

By  comparing  this  with  Eq.  (91),  we  see  that  p  has  the  meaning  of  the  operator 
of  the  angular  momentum  corresponding  to  the  internal  motions  of  a  molecule, 
and  it  is  known  to  be  called  the  operator  of  the  "vibrational  angular 
momentum." 

Equations  (91)  and  (92)  allow  us  to  rewrite  Eq.  (86)  as 
p(1)'  - 

and 

p.  .  '(*.»)  . 


(93) 

(93a) 


8.  Calculation  of  the  Operator  K’ 

When  calculating  the  RHS  of  Eq.  (85)  we  get  terms  both  quadratic  and 

% 

linear  in  Jtf-u  and  P.  The  former  ones  are  obtained  from  p'ep',  yielding,  in 
view  of  Eqs.  (93a),  (U9i),  (58),  and  (63), 

a, 8 

where  y  -  are  the  elements  of  the  matrix 

OLD 
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(95) 


r  -  G^r)'i^r)'. 

To  get  the  terms  of  K'  linear  in  tf-u  and  P  it  is  convenient  to  keep  the 
form  of  pv  '  in  Eq.  (85)  at  right  from  conana  as  given  in  Eq.  (93).  and  at  left 
from  conana  to  use 


p(l)'  - 


(96) 


Thereby  we  obtain  for  the  linear  terms,  in  view  of  Eqs.  (72),  (7h),  (58),  (77), 
and  (U9i), 

-  ^{i:ECjn^,G^r),AjrVn)^-Aa1A[r)^r^,(*f-p)  +  BjP] 


Ji 


-  lK.O(r)’l  (#-v)I 


-  -  4-  r(s«  ♦  s"K#-m)  -  %s.p. 

2  0  a  a  2j  J  J 


(97) 


(o)'  (a)' 

where  c'  '  is  Jth  row  of  Cv  •  . 

V 

Por  any  column  matrix  ■  (f^,  f4y,  fJz)  we  get,  in  view  of  Eqs.  (22), 
(2U),  and  (72), 

.<v(a) 

_  .(a) 

:«  «  *«  -  to*^  w 
Ji 


EIC<”)-G(‘-)Xr,*i‘,)  *3§-q  ■  E3pC(0)'o(r,')tjr,*5'J5c( 
ji  j  i 


-  SpG^'eX^A^I  (98) 

where  A*?*!  is  obtained  from  A[r*'  as  given  in  Eq.  (25)  by  substitution  f16 


for  5, 


Por  s  we  should  write 

Cl 


5  • 


(99) 
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where 


(100) 


is  oth  column  of  <j^'.  Eqs.  (30)  and  (100)  yield 

fi  ‘  ^ix’VV 

»  m  (R(i)K(r),-R(i)«(r)’  R(i)*(r),-R(iUr)’  »<*>>>’  „<*>>)'* 

“i  z  V  V  8zc  ‘Rx  gza  Rz  «xa  *  y  gxa  ~Rx  gya  h  (101) 

Furthennore, 


rX(r)A(r)f'» 

“i  Ai(f) 


where 


At) 


.(f) 

-I<f>  . 

■!<*> 

XX 

yx 

zx 

.(f) 

- 

xy 

yy 

zy 

.(f) 

-i<f> 

l(f) 

xz 

yz 

zz 

lut. 

)•  I(f) 

*  ER 

S  I 


(f) 


(102) 


(103) 


1(a#(Mt)  1 

In  what  follows  we  will  hoose  the  principal  inertia  axes  as  the  axes  of  the 
rotating  coordinate  system,  thus  obtaining 

0 

yza  z-ya  | 

0  1  (10U) 


At) 


I  g^' 
yza 

-i  gir)' 
z'Va 

,(r)' 

'za 

0 

I  g(r)‘ 
z6xa 

,(r) ' 
■ya 

0 

where  I  are  those  given  in  Eq.  (58a).  We  now  finally  get  froa  Eqs.  (99) 
01 

and  {10k) 


s'  ■  SpG*r)Vf) 


*08  *y6  '  nr*«Y  g6o  *@r  g«*  (ororv) 


(105) 
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(i)  (r)»  • 

For  small  enough  6r  ,  when  G  is  a  symmetric  matrix,  s  obviously 
on  the  form 


i*  *  (g(^Vr),-Y(r),g(r)’)(I  -I  )• 
a  'saf5  say  ’aa  y 


For  Sj  we  have 


(i)  (i)  (ik 

bJx  *bjy  ’V  > 


which  leads,  with  due  allowance  for  Eq.  (55c),  to 


W) 


0 

bii} 

-b<i) 

J* 

Jy 

h(i) 

■v 

0 

b<i) 

3* 

b(i) 

bJy 

„<i) 

bJx 

0 

where 


At)' 

*1(0 


H  (r)' 

WiU) 

t»l 


a<i) 

a(i> 

f 

aia 

*y 

it 

0 

a(l) 

ix 

,(i) 

iy 

-a(1> 

ta 

0 

Therefore,  in  the  case  under  consideration, 

N 


r^(r)A(r)'  -  r  >(r)A(r)f 
J*i  Ai(f)  r8Jif“iAi  AiU) 
1  1-1  1 


But,  in  view  of  Eq.  (68),  the  BBS  of  Eq.  (110)  can  be  rewritten  as 

J^U) 

where 


Hence,  we  finally  obtain 


tanes 

(106) 

(107) 

(108) 

(109) 

(no) 

(in) 

(112) 
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(113) 


where 


sj  ■  S*°M JV<«>  -  r«j,c« 

l-l 


N 
rg 

i-i 


(r ) '  ' 

Ct  -  SpG<r>  I(t). 


At  last,  as  follows  from  Eq.  (97) 


s"  «  -Ik ,G 
«  4  J 


a>(r) '  .^(r)' 

(r)  1^2 - SptcG  r  'i  -2— 


3*4 


»q 


But 


aS(r)'  .  aii^(r)* 


3q 


3q, 


Therefore,  in  view  of  Eq.  (112), 

3&(r)' 


**» 

aq 


(q) 


where 


and  Eq.  (llU)  with  due  allowance  for  Eq.  (95)  takes  on  the  form 


Combining  all  the  results  obtained,  we  get 

K’  •  i  *  2  Wj,y» 


a, 8 


-  T  £3ipi  "  ^  E(s«  +  8«)(wa"w«)- 
‘  j  J  J  2aa  a  00 


.(1). 


All  the  terms  of  Eq.  (119)  linear  in  M  -y  vanish  when  5r'  ;-K),  but 

a  a 

linear  in  Pj  do  not.  This  fact  can  cause  difficulties  in  practical 


(113a) 


(111*) 


(115) 


(116) 


(117) 


(118) 


(119) 

those 
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utilization  of  the  Hamiltonian  containing  K' ,  and  it  is  more  convenient  to 
replace  the  latter  by  an  operator  K",  auch  that  it  does  not  contain  the 
undesirable  terms  and,  at  the  same  time,  is  connected  to  K'  by  a  similarity 
transformation.  This  transformation  should  make  use  of  the  Jacobian  of  the 
initial  Cartesian  coordinates  with  respect  to  the  new  ones  ve  have  chosen, 
as  it  takes  place  when  passing  to  the  normal  coordinates. 


9.  Calculation  of  the  Jacobian 

The  Jacobian  of  the  initial  Cartesian  coordinates  with  respect  to  the 
new  ones  obviously  is  |0|  ,  where  I>  is  the  matrix  given  in  Eq.  (3M.  Prom 

Eq.  (3M  we  get 


ME.  0 

0 

* 

0  2f(Q)j*0(o) 

0  itsA(r),C{Q) 

T 

( 120 ) 


where 

I"  -X(r,Vr)\ 

Equation  (120)  yields,1^  in  view  of  Eqs.  (55c)  and  (76), 


(121) 


|8oD|  «  |D|2  «  eonst[C{n*|2|T|  [ I^-kTic | .  (122) 

(i)' 

Since  we  are  interested  only  in  the  factor  depending  on  r  ,  we  will  deal 
merely  with  | l"-rTic | .  To  simplify  this,  let  os  consider  also  the  inverse 
matrix 


D"1e2r1 


vT\  0  0 

0  c(ft),rfr(0)’  -c{n),r* 

0  G  ♦  rrr 


(123) 
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It  follows  from  DmD(D~1eD“1)  *  E^,  in  view  of  Eq.  (6l),  that 

i"-kT(c  =  r"1. 

Therefore,  in  view  of  Eqs.  (122),  (95),  and  (71b), 

| D (a 1 1 '  ( ;  |Dr1«|G(r),|. 


(121*) 


(125) 


1  I 

: 


10.  Calculation  of  K”  and  the  Hamiltonian 

In  this  section  we  will  need  some  results  concerning  canonically 

IT  18  *  * 

conjugate  operators.  *  That  term  is  used  for  any  operators  a  and  b 

obeying 

[ab]  ■  cE  (126) 

where  c  is  a  constant  and  E  is  the  unit  operator.  Operators  and  qj  are 
canonically  conjugate,  £  being  equal  to  -ifc.  Functions  of  canonically 
conjugate  operators  are  noncommutative ,  and  if  there  are  R  mutually 

A  A 

independent  pairs  of  canonically  conjugate  operators  a.  and  b.  for  which 
the  relations 

[aj\]  “  l*jV  *  *  0  (127) 

where  E  is  the  unit  operator,  are  valid,  then  for  the  operators  4  ■ 

A  AA  A  AA  AA  A 

^ « y $  •  •  •  }bj| )  cuid  f  *  f  (&^ $  •  •  •  ^a^  |b^  |  *  • » |b^ )  we  oht aln 


a-z* 

£  1  3b 1  3a  3b  3a* 


(128) 


where 


JL  iS  lli 

•  I  *i 

tl  3b*  3a* 


.(^...tjjj) 


-1 


•* 

3*f 


* ' 


,lR 


1* 

JJl 


"  1,  .1m  .  *i|» 

3b1A...3bN"  3a11...3a1|" 


(129) 
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When  applying  Eq.  (128)  to  the  calculation  of 


K"  ■  |G(r),r1/,2K'lG*r)’|1/2 

we  will  need  the  first  and  second  order  derivatives  of  |G^  l1^2  with 
respect  to  the  q^s .  We  have,  in  viev  of  Eq.  (112),  and  (113a) 


»is!.r),l^a ,  k  (G(r),(l/2SpI,is.l!l.  - 1  fG(r),«1/g- 


3q 


J 


and,  in  view  of  the  linearity  of  I*  in  the  firm's, 

£[0(r)’|l/2 


Gvr;  l1'*  ,  1  |f,(r),,l/2-  r(r)'  .  ^(r)*. 

—  t  |G  1  hi)^0  T(J) 


(130) 


(131) 


♦  \  |0<r),|1/2Sp0(r),ljt)0(r)'l(,3)  -  \  |o(r)’|1/2(|-cJlc3+cl3)  (132) 


where 

■s»0<rl'Iu)0<r,'I(j)'  (l3a,) 

In  Eq.  (130)  |o^r^  |1^2  is  a  function  of  the  6r'*^  only,  whereas  the 
U  are  functions  of  both  the  P,  and  fir^.  But,  as  it  is  seen  from  Eq.  (128), 

<*  J 

we  will  not  need  the  derivatives  of  u  with  respect  to  q,,  since  |G^  | 
does  not  depend  on  Pj.  On  the  strength  of  Eqs.  (131)  and  (132),  the  trans¬ 
formation  given  in  Eq.  (130)  merely  adds  to  K'  sone  linear  terms  and 

functions.  In  particular,  the  linear  terms  resulting  frosi  j  Z  *re. 

J  ,t 

in  view  of  Iqs.  (128),  (113),  snd  (131), 

■  f  ]-iTi  ■  ,133> 

These  obviously  cancel  the  corresponding  linear  terms  in  K*.  The  term 
quadratic  in  Pj  yields  also,  in  viev  of  lq.  (132), 
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I 

l 


vi  *  ~  r r  gij(2  Vj +  cij}- 
From  the  term  linear  in  P.  we  get,  in  view  of  Eq.  (131), 

.  *  tf 

VI  *  r  ;»JCJ  -  T  E  *t3ctc y 


( i  :;'i ) 


(1 35) 


Ve  get  linear  terms,  which  are  due  to  the  term  quadratic  in  #  -w  ,  from 

a  a 

ill 


x  lW.c3  w. 

n  a  J  J 


0,8 

In  view  of  Eqs.  (92),  (7!*a),  and  (128), 


(136) 


- 5p/'f  ■  -  'W'n 


“  ~  ♦  *«*(*<.”>*  )J  “ 


Jo'  p  “S'  J8  o  -o 


i  J*  te 


(137) 


where 


(U)  -  fill 
»«t  * 


(138) 


Therefore,  Eq.  (136)  yields 


o8  *  CJ*JB  *** 


2  1  Yo8?° JK JB^a”wa ^  +  V  1  Ya8  1  " 
o.B  J  a,B  J,i 

ill 


-f  js:w  +  v2 

a 


(139) 


where 


s'"  «  £Ec  ic  y 
a  J  JB'Bo 


(lUo) 


and 


42  /a\ 

V2-TE  "jl'taW 

o.B  i,J 


(1  kl) 
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From  the  same  quadratic  term  we  also  obtain,  in  view  of  £qs.  (126), 

(132),  and  (137), 

*2  . 

V3  *  "  T  1  Yo0  1  <*}«*«  +  *cta'eJB^2  ctcJ  +  cij} 
a,B  J,t 

*  "  T  Z  Z  Ya0,cJa,ei8*2  °iCJ  +  C4J)*  (lt2) 

a.B  j,t 


At  last,  terms  linear  in  M  -v  yield 

a  a 


\  -  -  V  IK  * 


(1U3) 


Combining  all  the  above  result*,  we  can  finally  write  the  rotation- 
vibration  Hamiltonian  H"  *  K"  ♦  V  as 

"•■I1  VV'.,aW  *  1 1  «jiY« 

a.B  J,* 


where 


and 


-  i£ 


2  Ia«(W  +  V’  +  V 

0 


s  *  s'  ♦  s'*  ♦  s'" 
a  a  n  a 


V'  ■  VJ  ♦  VJ  ♦  v2  ♦  v3  ♦  vh 


(lkk) 


(H»5) 


.U) 


■*•??  .:/5rv- 


*  r  *2  gJt  “  1  *.Ja*MY«^cJc*  "  1  *eJi  *  2  1  KJaKtBYo0^cJt^  (1^ 

.)»!  a, I  J*1  ®,0 
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Here  s^,  s£,  and  s£  are  as  given  in  Eqs.  (105),  (118),  and  (lfco),  respective¬ 
ly.  The  rest  of  notation  renalns  the  same  as  before. 

11.  Discussion 

The  coefficients  Ya8  and  sa  in  Eq.  (lUU),  as  well  as  V*  in  it,  can  be 
expressed  explicitly  as  functions  of  Cartesian  displacements  of  atoms  in  the 
rotating  coordinate  system  and  do  not  depend  on  the  choice  of  internal 
coordinates.  Therefore,  they  can  be  calculated  in  a  system  of  internal 
coordinates  the  most  appropriate  for  every  particular  case,  irrespective  of 
what  a  system  of  internal  coordinates  could  be  chosen  for  the  computation  of 
molecular  vibrations. 

lo  restrictions  in  the  magnitudes  of  the  atomic  displacements  have  been 
made  in  deriving  the  new  Hamiltonian;  therefore,  it  is  exact.  In  contrast 
with  this  Hamiltonian,  the  conventional  one  is  valid  only  for  small  enough 
atomic  displacements  due  to  the  fact  that  they  are  represented  as  linear 
combinations  of  normal  coordinates. 

The  nev  Hamiltonian  can  be  easily  represented  in  the  form  of  a 
differential  operator,  by  making  use  of  Eqs.  (91),  (92),  (71a),  and  (5). 

This  fact  releases  one  from  the  necessity  of  utilizing  the  only  computa¬ 
tional  method  used  by  now,  namely,  the  perturbation  one,  and  allows  one  to 
hope  that  some  other  methods  of  approximate  calculation  could  be  used  more 
successfully.  Even  if  one  chooses  the  perturbation  method,  the  new  Hamil¬ 
tonian  seems  to  be  simpler  for  applications,  since  It  has  an  explicit  form 
and  all  its  coefficients  can  be  expressed  explicitly. 

In  application  of  the  new  rotation-vibration  Hamiltonian  to  the 
particular  case  of  linear  molecules,  if  the  molecular  axis  is  chosen  as  the 
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z-axis  of  the  rotating  coordinate  system,  then  x  *  0*  the  and  R^‘^  are 

zeroa  in  Eq.  (30)  and  everywhere  el3e,  because  of  which  I  and  I'  as  well  as 
Gv  7 ,  Gv  ,  and  T  are  2  by  2  diagonal  matrices  with  equal  diagonal  elements 
in  each  of  them.  The  Hamiltonian  as  given  in  Eq.  (lUU )  is  thereby  signifi¬ 
cantly  simplified. 

I  will  emphasize,  too,  that  in  deriving  the  new  Hamiltonian  the  only 

supposition  that  has  actually  been  made  is  the  existence  of  an  equilibrium 

19 

configuration  for  a  molecule.  No  Eckert' s  conditions  proved  to  be  necessary 
anywhere  in  the  deduction. 

In  the  process  of  the'  derivation,  the  equations  of  fU  proved  to  be  very 
helpful.  The  method  of  calculation  of  the  Jacobian  needed  (f9)  as  well  as 
the  method  used  for  the  similarity  transformation  of  X'  (§10)  may  be  of 
particular  interest  too. 

The  following  advantages  of  the  new  rotation-vibration  Hamiltonian  over 
the  conventional  one  are  obvious:  (l)  the  exactness  of  the  former,  in 
contrast  with  the  approximate  nature  of  the  latter;  (2)  the  explicit  form  of 
the  former,  in  contrast  with  the  tangled  fora  of  the  latter;  (3)  the  indepen¬ 
dence  of  all  the  terms  of  the  former  but  V  of  the  force  coefficients,  whereas 
all  the  terms  of  the  latter  except  V  are  implicit  functions  of  those;  (M  the 
possibility  of  applying  approximate  calculation  methods  other  than  the  pertur¬ 
bation  one  to  the  former,  whereas  that  is  impossible  to  do  for  the  latter. 

Certainly,  scam  other  advantages  of  the  new  Hamiltonian  over  the  conven¬ 
tional  one  would  be  revealed  in  practical  application  of  the  former  to  the 
computation  of  various  particular  rotation-vibration  spectra.  The  new 
rotation-vibration  Hamiltonian  can  prove  to  be  of  special  importance  for 
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computation  of  molecular  spectra  in  the  case  of  large  rotation  or  vibration 
quantum  numbers,  in  particular,  for  computation  of  spectra  of  molecules  near 
their  dissociative  limits.  But  I  do  not  foresee  cases  vhen  the  conventional 
rotation-vibration  Hamiltonian  could  have  any  advantage  over  the  new  one. 
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